Solution of mass/gap equations for strong velocity anisotropy in the QED3 theory of 

underdoped cuprates 
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The low-lying excitations at the nodes of the d-wave gap in the normal state for underdoped 
cuprates, close to the superconducting phase transition, may be described by an effective QED3 
theory. There are three characteristic velocities: vp, va and c. For ua = = c, the model reduces 
to A''-flavour QED3. Here, in the isotropic limit, for A'' < A'^^, a critical number of fermions, a 
dynamical mass is generated which corresponds to the formation of a spin density wave. We study 
the effects of strong velocity anisotropy {vf 7^ va) on dynamical mass generation. Solutions are 
given for the dynamical mass at both N = 2 and A'' = 100, and so we show that the critical number 
of fermions Nc > 100. However, we argue that Nc — N^, when we go beyond the approximations 
used to derive our mass gap equations. We expect, though, that our solution for A'^ = 2 is roughly 
correct, at low momentum, for large enough anisotropies. Implications and possible extensions to 
this work are also discussed. 



I. INTRODUCTION 

Based on two simple assumptions, recently, it has been 
possible to derive an effective theory jl] , |^ that describes 
the low energy excitations of the cuprates in the under- 
doped regime. Firstly, that there is gap of d-wave sym- 
metry, to insure that there are gapless quasiparticle ex- 
citations at four points on the Fermi surface as well as 
a gap that changes sign. Experimentally, it has been 
reasonably well established that these materials possess 
d-wave symmetry. The results of in particular, show 
a linear dependence in the penetration depth with re- 
spect to temperature, suggesting that there are, indeed, 
nodes in superconducting gap. In Q| there is further 
evidence to suggest that this assumption is correct; an 
observed change of sign with direction, in the supercon- 
ducting gap, which is consistent with d-wave symmetry. 
In the underdoped regime, the cuprates can be consid- 
ered as two dimensional Therefore, as a second as- 
sumption, one could reasonably expect that such materi- 
als undergo a Kosterlitz-Thouless type phase transition 
1^, in which vortex-antivortex pairs destroy the super- 
conductivity. Already, there is some justification for this 
second assumption, through the experimental evidence 
of and [D , that supports the idea of free vortices near 
and above the superconducting critical temperature. 

The effective theory of Q , [|| is an anisotropic version 
of 2-1-1 dimensional quantum electrodynamics (QED^), 
with three characteristic velocities; vp the Fermi velocity, 
VA a velocity associated with the steepness of the gap at 
the 'nodes', and c the propagation velocity of the gauge 
field. Here, the isotropic limit is realized when vp = va = 
c. The theory consists of two {N — 2) species of Dirac 
fermion that correspond to neutral spinous near the two 
pairs of diagonally opposed nodes in the d-wave state, 
as well as the gauge field that describes the interaction 
of the spinous with vortices. Spinous can be connected 
back to the real quasiparticle excitations in the cuprates 
through the singular gauge transformation discussed in 



In jl] it was shown that at T = there is inherent 
instabi ity towards an antiferromagnetic state when one 
has unbound vortex loops in such a model. In QED3 this 
instability corresponds to a " chiral" instability in which a 
mass is generated for fermions, dynamically, though their 
interaction with the gauge field. In the isotropic limit, 
such an instability is present when the number of Dirac 
components N < Nc, with iVc w 3 This number 
seems to be fairly robust, so as to stay roughly the same 
when more sophisticated approximations are employed 
Ip^ , ||ll[ . In the context of d-wave superconductivity, this 
dynamically generated mass is in fact a staggered poten- 
tial felt by the original electrons, i.e. weak spin density- 
wave (SDW) order ||l|. One should point out that such 
a spontaneously generated SDW automatically confines 
the neutral spin- 1/2 excitations (spinous) of the super- 
conductor pa , so ruling out spin charge separation when 
unbound vortices are present. 

It is important to understand how velocity anisotropy 
in QED3 affects the generation of a mass/gap. This 
problem is not only relevant to SDW generation in the 
cuprates, it also appears in the other effective Dirac the- 
ories describing condensed matter systems jl^, which 
are not required to be Lorentz invariant. In the 
effect of weak anisotropy was examined in the large-N 
approximation of |9|, with two main results. Firstly, that 
weak anisotropy becomes irrelevant when the mass/gap, 
S — !■ 0; this was found to be in agreement with a similar 
study in 'chirally' symmetric state Secondly, it was 
argued that Nc is the same in the weakly anisotropic case 
as it is in the isotropic case, if, indeed, Nc can be consid- 
ered universal. For the cuprates, however, the anisotropy 
is quite large vp/vA ^ 10, p6) . So indeed one might ask 
what the effects of strong anisotropy might be. 

We shall set out to answer such a question by first 
solving the gap equation in the strong anisotropic limit. 
Here, we employ the same approximation scheme as [ p^ , 
but now we focus on the strong anisotropic limit. For our 
solution to the mass gap equation we find that Nc > 100. 
Secondly, we find for all our solutions that the mass / gap 
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intially increases with momentum, this effect becomes 
more pronounced as the anisotropy is increased. We 
then go a step further and calculate the wavefunction 
renormalizations for each of the directions in momentum 
space. We find that v4i « 1, where Ai is the wavefunc- 
tion renormalization in the direction most important for 
dynamical symmetry breaking in the strong anisotropic 
limit. This tells us that the approximation scheme we 
employ is likely to be reasonably accurate, provided that 
N < Nc, for strong anisotropics. We also calculate the 
renormalized anisotropy as a function of momentum, and 
discuss its effect on the validity of our approximation for 
strong anisotropy. Numerically, we find that Nc > 100, 
this may suggest either that Nc is very large or infi- 
nite. If Nc — > oo, we expect that this is an artifact the 
strong anisotropic limit, and that Nc — > oo only when 
the anisotropy tends to infinity. Furthermore, provided 
that Nc is finite, we do not suspect that N'^, the critical 
number of flavours calculated when the effects of renor- 
malized anisotropy are considered self consistently in the 
Swinger-Dyson equations, is as large as Nc- Therefore, 
we would suggest that, for strong anisotropics, at some N 
very close to N^ our approximation scheme breaks down. 
Based on two assumptions, Nc being finite and univer- 
sal, we argue that A^^ = iV°, the value calculated in the 
isotropic limit. 

The layout of the paper is as follows. In Section II we 
derive the mass gap equation in the strong anisotropic 
limit and proceed to calculate an approximate solution 
using Fade approximants. In Section III, we derive equa- 
tions for wavefunction renormalizations and the renor- 
malized anisotropy in the strong anisotropic limit. Here, 
we calculate the effective values of the wavefunction 
renormalizations as well as the renormalized anisotropy. 
Finally, we present the argument why it could be the case 
Nc = N^, its value in the isotropic hmit. In Section IV, 
the disscussion, we discuss the limitations, and possible 
implications of our results. 



II. DYNAMICAL MASS GENERATION IN THE 
LIMIT OF STRONG ANISOTROPY 

The effective low-energy theory for the d-wave quasi- 
particles (i. e. spinons) at T = coupled to the topolog- 
ical defects (vortex loops |jl^ ) in the superconducting 
phase can be written as £ = £i + £2 + with |l|: 



N/2 



70(^0 - ieao) + -j=-f2{d2 - «ea2) 



£a 



(5VA7i(92 — ieai) 
(V X a)2 1 



2e 



5(a) (1) 



and £2 is the same as £1, except that l/-\/A ^ a/A. Here, 
5 = ^JvfvaIc and A = vp/v^- The 7-matrices satisfy 



the Clifford algebra {7^,7;^} = —S/^^. The gauge field 
a represents the singular superconducting phase flucta- 
tions induced by the vortex loops, and VE" fields represent 
the spinons; the neutral spin-1/2 excitations one can de- 
fine near the nodes in the dSC. The coupling constant 
e ("charge") is proportional to the dual order parame- 
ter that signals the appearance of infinitely large vortex 
loops, i.e. the destruction of the phase coherence in the 
dSC fill, @- We have assumed N/2 identical copies of 
each type of spinon field; the value of physical interest for 
a single layered superconductor being N=2. Here G{sl) 
is a generalized gauge fixing term. We shall work in the 
landau gauge, as in Q , so one chooses G — (V.a)^ and 
e = 0. In the limit S ^ 1, A ^ 1 the model simply 
reduces to the much studied iV-flavour QED3. 

When studying dynamical mass generation in the 
isotropic limit, one may simply approximate the effect 
of Maxwell term simply by replacing (V x a)^ by a u.v 
cutoff on the momentum integrals A « Ne"^. Essentially 
this can be done, because the effect of the Maxwell term 
is to cause E to fall off rapidly at p = A, the momenta 
at which maxwell term becomes significant, so effectively 
cutting of the integral. But, considerable care must be 
taken in the strong anisotropic limit; the Maxwell term 
can no longer be neglected in such a way. We shall show 
that it is needed to regularize a i.r divergence present in 
the ko integral if one uses solely the polarization tensor 
in the strong anisotropic limit. Since, now, the Maxwell 
term can no longer be neglected it remains unclear what 
the role of 5 is in the large A limit. This is because any 
non-trivial dependence on S can no longer simply thought 
of as an artifact of gauge fixing , for scaling 6 out of 
the action rescales both the gauge fixing term and the 
Maxwell term, which we are forced to retain. However, 
now for the sake of simplicity, in Eqn. 1 we shall assume 
6 = 1, and the problem oi S ^ 1 will be left to a future 
publication. 

Now, we shall calculate the mass gap for the type 
1 spinons (those a appearing in £1) in the strong 
anisotropic limit. As discussed in jl^ , due to the pi —> p2 
symmetry of the theory, it is very easy to relate this to 
the mass gap for the type 2 spinons. From Q, we recall 
what the form of the full spinon propagator should take: 



(2) 



SR{k) = {j:ik) + A^{k)-/^k^) 



Here, the repeated index, as usual, implies a sum over 
indices. S is the dynamically generated mass-gap. 
is refered to as wave-function renormalization in k^ di- 
rection. The starting point in our calculation of S is the 
form of the gap equation given in p4| with 5=1 



(Pk I](fc) 

(2^)3 



where 



U{k) 



fcg 



Afcf 



Si(fc)^ 



(3) 
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(4) 



with q — k — p,k — (fco, ki,k2), and a = Ne^. Using this 
as weh as the results given in for the polarization 
tensor and the full gauge field propagator to one loop 
order, with massless spinous, one is able to extract the 
following limiting form for a gap equation in the Landau 
gauge. 



16\/A [ (fk E{k) 



N 



(2^)3 \kf + S(fc) 



n{q), 



(5) 



where 



f 162glg2g^ , 



k|2(go +g2 



16(|g2|(gg+gf) + |gl|(gg+g|)) ^ 



(6) 



It is not hard to see that if the terms arising from the 
maxwell term (terms multiplied by l/aV\ and in 
Eqn.6 ) were neglected then one would, indeed, have a 
divergence in the fco integral at po. Next, we proceed by 
shifting the ko and fc2 integrals so that H is independent 
of Po and p2- Then, it is easy to show that there is a 
solution to Eqn.5 of the form T,{p) = M(pi)aA/16 where 
Pi = Pi/Cc^a/A). On rescaling the momentum integrals 
we may write 



M(pi) = 



VXn^N Jo 



dki 



M{ki 



kl 



M{ki 



7^(91), (7) 



where Z^(gi) = / dko J dk2'H{ko,qi,k2), and now gi = 
kl — pi- We first perform the fco integration; then when 
we come to performing the fc2 integration, we see that 
some of the terms contain logarithmic u.v divergences in 
fc2- To cure these divergences we first subtract out the 
singular part of U{qi). We are able to show that the 
singular part of U (gi ) takes the following form 



Z^(gi) 



Sing 



G(gi) 



dk2 



1 



l + fca 



where 



G(gi 



0.375|gi| +0.2571 
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(8) 



(9) 



These divergences should be considered an artifact of ne- 
glecting the fc| term in U{k). Therefore, to render our 
mass/gap equation completely finite we should then re- 
place U{qi) sing (kj + Af(fci)2)-i with 



U{qi)cor = G(gi) 



dk-2 



1 



1 



l + fc2 fcf + fci/A2+M(fc-i)' 



(10) 




FIG. 1: Graph showing M calculated for = 2, A = 
10, 15, 20, 100, as a function of pi, the rescaled momenta. The 
curve with the smallest magnitude is for A = 100 and the 
largest curve is for A = 10 , the curves decreasing in size with 
A 



On evaluating this integral for large A we are then able 
to recast our gap equation in a final form 



M(pi 



Af(fci) 



VXtt^nJo fc? + M(fci)2 



J^(fcl,gl) = [U{ql)f^^ 



(fc2+M(fci)2)l/2 



•^(fcl,gl), 

log(A)G(gi)+G(gi) 



-l0g(fc2+M(fci)2) 



(11) 



where U{qi)fin ^U{qi) -U{qi)sing- The ^2 integral for 
^(gi) fin needs to be evaluated numerically, we find that 
^(gi) fin is reasonably well approximated by 



U{q)fr 



(.01357 + 62.84|g| - 39.04g2 - 5.476|g|3) 
(.3109 + 88.55|g| + 78.61g2 + 3.486|g|' 



(12) 

As a further approximation we also replace M(fci) with 
M (0) in all the dominators as well as the logarithm, an 
approximation valid for M < I. Then, we calculate an 
approximate solution for the mass gap equation using an 
approximate solution of the form 



M{p) = 



mi + 7712 + 7773p2 
1 + 7774 + 7775J32 



(13) 



then by calculating M(0),M(1/3),M (2/3), M(l) and 
M(5) for trial values of rrii {i = 1..5) one calculates new 
values of rrii. We then iterate the process until the M (p) 
inserted into R.H.S is the very close to the M{p) calcu- 
lated for L.H.S in Eqn.ll 

In Fig.l we show the calculated values of the mass/gap 
for = 2 for A = 10, 15, 20, 100. As A is increased the 
solutions are seen to fall in magnitude. One can see that 
all the solutions seem to possess three qualitative fea- 
tures. The first is that for small momenta we see that 
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FIG. 2: Graph showing M calculated for N = 100, A = 
10, 100, as a function of pi, the rescaled momenta. The curve 
with the smaller magnitude is for A = 100 and the larger curve 
is for A = 10. 
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FIG. 3: Graph showing A\ calculated for A'' = 2, A = 
10, 15, 20, 100 in units of \/A, as a function of pi , the rescaled 
momenta. The smallest curve is for A = 10 and the largest 
curve is for A = 100 , the curves increasing in size with A 



the mass/gap increases with rescaled momenta. Then a 
maximum in M(j)\) is reached and then the solutions de- 
crease with momenta. Finally, we see that the solutions 
tend off to a constant for large momenta. However, as we 
discuss later in section IV, this last feature may well be 
an artifact of our approximation. As A is increased the 
first and third feature become more pronounced, whereas 
the solutions seem to drop off less for momenta above the 
maximum in M . 

In Fig. 2 we show the calculated values of the mass/gap 
for = 100 for A = 10, 100. Again, we see that as A is in- 
creased the solutions are seen to fall in magnitude. How- 
ever, we see that no-longer a maximum in M can be seen; 
these solutions seem to be monotonically increasing with 
momentum towards some constant value. Since a solu- 
tion exists for A^ = 100 we may conclude that > 100 
for Eqn. 11, where only bare anisotropics are considered. 



III. THE WAVEFUNCTION 
RENORMALIZATIONS AND RENORMALIZED 
ANISOTROPY 

From [|l4[ we have the following expressions for the A,, 



the (p) depend only on pi and 



A,{p)^Xjl~ — 



where 



T^(fc,g) = koD^^oiq) + \k,D^M + tl^^M, (15) 

Here, Aq = 1, Ai = a/A and A2 = I/a/A- From these 
expressions one is able to show in the large A limit that 



AM^Vxil 



Aoipi) = VXA2{pi) = 1 



16\/A f (Pk k{H{q) 
Npi J (2^)3 Xkf + j:{ki 



16VA f d^k 



n{q) 



N J (27r)3 Afc2 + E(/!ci)2 



(16) 



where Ti{ki) = I](0, fci, 0). Then it is easy to see, that to 
calculate the , one goes through the same analysis dis- 
cussed in the previous section. We arrive at the following 
expressions from Eqn. 16 



^o(pi) = l 

Vx 



dki- 



T{ki,qi) 



VXirm Jo kl + M{k 



dkl ^'^tn^l (17) 



VXTT^NpiJ-^ kl + M{ki) 



In Fig.3 we calculate Ai for A^ = 2, A = 10, 15, 20, 100; 
also the renormalized anisotropy A^. = yfXAiipi) / A^ipi) , 
again for A^ = 2, A = 10, 15, 20, 100 in Fig.4. In exam- 
ining Fig.3, the first thing to notice is that Ai does not 
depend much on pi and is close to 1. This means that 
the effect of including vertex corrections, as well as wave 
function renormalization for the fermion propagator, in 
the self energy might be very small. This is assuming the 
vertex ansatz F^ = A^"f^ (or some other ansatz depend- 
ing only on A^) is a good approximation, where F^ is the 
full vertex function. This ansatz is simply a generaliza- 
tion of the ansatz F^ = A'j^^ used in the literature when 
considering the isotropic limit, jll] ||T^. This assumes 
that S can be neglected in vertex corrections. 

In Fig.4 we see that the renormalized anisotropy is ef- 
fectively a half of its bare value, and changes little with 
momenta. Although this does effect the equation for the 
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self energy, which depends only on Ai and S in strong 
anisotropic limit, this may well effect the polarization 
tensor, though by how much is not clear. This is simply 
due to the fact that all elements of the polarization ten- 
sor, with the exception of Ili^2i contribute in the strong 
anisotropic limit. This implies that all the contribute. 
However, such an effect on the polarization may lead to 
a small change in E. More importantly. A,- will effect at 
what values of A our large momentum approximation is 
valid. The strong anisotropic limit should be considered 
to be a good approximation above some specified value 
A = As, so far undetermined. The fact that the renor- 
malized anisotropy is approximately half the bare value 
tells us that As should be effectively doubled, when the 
renormalized anisotropy is included in the self energy. 

Now we present an argument that Nc{X) = Nc{X — 1). 
There are two assumptions that we shall make. The first 
is that iVc(A) is universal; that Nc{X) depends only on 
Xr{0) in the limit M(0) 0. The second is that for 
finite A there exists a 'bare' (calculated with only the bare 
anisotropy A) finite Nc, which monotonically increases so 
that Nc ^ OO as A ^ cxD. Such a finite Nc probably arises 
from terms neglected in the strong anisotropic limit, as 
our result Nc > 100 suggests that Nc ^ oo for Eqn.ll. 
At such an Nc, M(0) = 0. If we are able show that for 
M(0) — > that Xr flows towards zero, then, the effect 
of a faffing A is to reduce Nc- We know from Q, 
that Xr{0) flows towards 1 for weak anisotropics, so if it 
does so for strong anisotropics then it is probably safe to 
assume that A is irrelevant for all values. So, assuming 
that these assumptions are correct, then indeed it must 
be the case that Nc{X) — Nc{X = 1). We are able show 
that Ar(0) is irrelevant in the strong anisotropic limit. 

To do this wc must examine both expressions in Eqn.l7 
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FIG. 4: Graph showing calculated for TV = 2, A = 
10, 15, 20, 100, in units of A, as a function of pi, the rescaled 
momenta. The curve with the smallest magnitude for large 
rescaled momentum is for A = 100 and the largest curve is 
for A = 10 , the curves decreasing in size with A for large 
momenta 



in the hmit M(0) 
we notice that 



0, for pi = 0. In the limit pi — > 



Va 



= 1 



^/Xtt^N 



ki 



dki- 



dqi 



Pl=0 



ki+M{klY 



(18) 



Because of the singular i-r behaviour of both expressions 
in the limit M(0) — > 0, we require only the small qi be- 
haviour of both !F{ki, qi) and its derivative with respect 
to qi. By extracting the small qi behaviour of J-{ki,qi) 
we are able to show that the dominant singular behaviour 
in the limit M(0) ^ is 



w \/A ( 1 
1 



7r2yAM(0)3/2 ) ' 
4 



(19) 



7r2VAM(0)3/2 

From these two expressions we may assume for A >> 1 



A„ ~ A 



7r2M(0)3/2 



(20) 



where A' = A — 1 and AJ, = A,.(0) — 1 (this is certainly true 
to 0{l/N)). Then Eqn.20 suggests that A; ~ M(0)3 , so 
we conclude X' decreases towards zero. 



IV. DISCUSSION 

In Section II we calculated the mass gap equation for 
TV = 2 and A = 10, 15, 20, 100, as wefl as iV = 100 and 
A — 10, 100. An important feature, pointed out in Sec- 
tion III for these solutions, is that for large momentum 
m{pi) tends to a constant. This feature is likely to be a 
artifact of our approximation pg| ] , and in reality we might 
expect our solutions to fall off slowly at large rescaled 
momenta pi for large, but finite A. However, we expect 
that our solutions are reasonable approximations at small 
rescaled momenta. What is reassuring in regards to this 
conclusion, is that we find that the wave function renor- 
malization Ai calculated in Section HI is close to one. 
As discussed there, this we would take as a sign that 
the approximation scheme we employ to derive Eqn.5 , 
our intial mass/gap equation, is reasonably accurate in 
the regime of strong anisotropy. So our conclusion would 
be that the feature of the mass gap intially increasing 
with momentum pi should certainly be true for QED3 
for large A. 

We should point out that the situation could be com- 
pletely different in the cuprates. Since the mass/gap is 
no longer cutoff at large momenta by the Maxwell term 
in our approximation, one would expect that some higher 
order term like the operator (V^ x a)2/e2 would cut the 
integral off at some scale pi ~ 1/{NX), thereby leading to 
significantly different behaviour. Such terms are present 
in the full effective theory |2^,[^ for the cuprates, once 
the vortex order parameter $ and auxiliary gauge field. 
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through which it couples to a, have been integrated out 
in the non-superconducting phase. 

Our most significant result seems to be that we can 
provide evidence for Nc{X) = N^, provided that our as- 
sumptions are correct. Namely, that Nc is universal and 
there exists a finite Nc, at finite A for E calculated with 
only the bare anisotropy in the mass gap equation. Even 
if the first assumption that we make is not correct for 
large anisotropies, we should take the fact that de- 
creases as M(0) 0, as an indication that Nc{X) changes 
little with A. The result that Nc > 100 suggests that it 
could be that iVc — > oo, but we should take this as Nc for 
the limiting case where X —^ oo. Detecting a finite Nc, or 
not, for bare anisotropy is likely to go beyond our strong 
anisotropic approximation, so requiring a more careful 
analysis of Eqn.5. However, it would seem strange that 
the second assumption, a finite Nc, did not hold. This 
would imply a Ac above which Nc — > oo. We think this 
is unlikely. The argument is unaffected by the inclusion 
of higher order operators, for example (V^ x a)^/e^, as it 
relies on the form of the small qi dependence of J-{ki, qi), 
where such terms can be safely neglected. The value that 
Nc takes for strong anisotropies has profound implica- 
tions in the model discussed in ||2l| . If Nc was to exceed 
10 for anisotropies where A ^ 10, might it be possible in 
some way to have coexistence of both antiferromagnitism 
and superconductivity. If our assumptions hold then the 
argument given in |21| , that there one phase transition, 
must hold for strong anisotropies. The conclusions, here 
again, suggest that for a single layered cuprate {N = 2) 



that there is no intermediate phase , between the 
SDW and the superconducting phase at T = 0. 

Clearly, an interesting problem for future considera- 
tion is the effect of higher order operators in the gauge 
field of 0(fc*) or higher on the solution of the mass gap 
equation. Even, in the isotropic limit the quantitative 
effect of these terms might be significant, even though 
such terms are unlikely to change qualitative behaviour 
of the mass/gap. For instance, such terms might affect 
the value of the momentum at which the mass/gap is 
effectively cutoff at = 2, but such a value will still 
be proportional to e^. As discussed before, the effect of 
such terms for large anisotropy is likely to become more 
pronounced, possibly leading to major qualitative differ- 
ences. Another interesting problem is what effect 5^1 
might have on the solutions to the mass gap equation for 
large A, this again is a problem we would like to address. 
There also remains the tricky issues of the universality of 
Nc and the form of the exact solution of Eqn.3, valid for 
all A, but these are much harder problems. 
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